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We study and quantify gravitational redshift by means of relativistic ray tracing simulations of emission lines. The emit- 
ter model is based on thin, Keplerian rotating rings in the equatorial plane of a rotating black hole. Emission lines are 
characterised by a generalized fully relativistic Doppler factor or redshift associated with the line core. Two modes of 
gravitational redshift, shift and distortion, become stronger with the emitting region closer to the Kerr black hole. Shifts of 
the line cores reveal an effect at levels of 0.0015 to 60 % at gravitational radii ranging from 10 5 to 2. The corresponding 
Doppler factors range from 0.999985 to 0.4048. Line shape distortion by strong gravity, i.e. very skewed and asymmetric 
lines occur at radii smaller than roughly ten gravitational radii. Gravitational redshift decreases with distance to the black 
hole but remains finite due to the asymptotical flatness of black hole space-time. The onset of gravitational redshift can be 
tested observationally with sufficient spectral resolution. Assuming a resolving power of ~ 100000, yielding a resolution 
of ~ 0.1 A for optical and near-infrared broad emission lines, the gravitational redshift can be probed out to approxi- 
mately 75000 gravitational radii. In general, gravitational redshift is an indicator of black hole mass and spin as well as 
for the inclination angle of the emitter, e.g. an accretion disk. We suggest to do multi-wavelength observations because all 
redshifted features should point towards the same central mass. 
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1 Introduction 

General Relativity (GR) teaches us that electromagnetic ra- 
diation is influenced by masses in both, photon energy and 
spectral flux. This gravitational redshift effect is remarkably 
strong for compact masses such as black holes and responsi- 
ble for their blackness. It is well established that black hole 
candidates in the universe hide in X-ray binaries (XRBs; 
Bolton 1972; van der Klis 1994, 2000) and active galac- 
tic nuclei (AGN; Lynden-Bell 1969, 1971; Rees 1984; Net- 
zer 2003). The existence of intermediate-mass black holes 
(IMBHs) in globular clusters (Baumgardt et al. 2003; Geb- 
hardt, Rich & Ho 2005) and ultra-luminous X-ray sources 
(Fabbiano 1989; King et al. 2001; Pakull, Grise & Motch 
2005; Miller 2005) is still under debate. Of course, gravi- 
tational redshift plays a major role in all these black hole 
mass families. According to the AGN standard model, there 
exists the so-called broad-line region (BLR) that is located 
between ~ 10 3 to ~ 10 5 r g away from the center. Here the 
gravitational radius is defined as r g = GM /c 2 with New- 
ton's constant G, vacuum speed of light c and black hole 
mass M (we set G = c = 1 throughout the paper). Doppler 
broadened optical emission lines with widths of typically 
10 3 -10 4 kms- 1 form in the BLR (Woltjer 1959). It has 
been first speculated by Netzer (1977) that gravitational red- 
shift may influence optical lines causing line asymmetries. 



Meanwhile, this has been reported in a number of sources, 
e.g. Peterson et al. (1985), Zheng & Sulentic (1990), and 
Kollatschny (2003). The gravitationally redshifted iron K 
fluorescence lines are very prominent X-ray features and 
have been reported for AGN and XRBs (Fabian et al. 1989; 
Tanaka et al. 1995; Martocchia et al. 2002; Miller et al. 
2004). First calculations of such lines around a maximally 
rotating black hole were performed by Laor (1991). A re- 
cent study by Nandra (2006) supports the idea that broad 
relativistic iron K lines are blended by non-relativistic nar- 
row lines originating most probably from the cold dust torus. 
The Schwarzschild factor allows for immediately estimat- 
ing gravitational redshift for static emitters. However, it is 
poorly understood how this effect evolves with distance to 
the black hole if a reasonable emitter kinematics is assumed. 
In this paper, we study the gravitational redshift over a large 
range of distances from the central black hole. We quantify 
the relativistic gravitational redshift on emission lines until 
GR fades beyond the current observable limit. The investi- 
gation is carried out in a very general form by discussing 
the observed line profile as a function of the generalized 
GR Doppler factor (g-factor) for Kerr black holes. Details 
of this work are presented in Miiller & Wold (2006). Previ- 
ous work on the that topic was done by Cunningham (1975) 
and Corbin (1997). 
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2 Method and model 

We compute emission lines with the Kerr Black Hole Ray 
Tracer (KBHRT) that maps emitting points in the equatorial 
plane of the black hole to points on the observer's screen. In 
a further step, we calculate the spectral line fluxes by numer- 
ical integration. All relativistic effects such as gravitational 
redshift, beaming and lensing are included, but higher order 
images are not considered. The solver has been presented in 
Miiller & Camenzind (2004, MC04 hereafter). 

We analyse the simulated line profiles by considering 
the ^-factor which generally reads 

.9 = fobs/^cm = Aem/Aobs = Z — ; , (1) 

1 + 2 

where v and A denote frequency and wavelength, respec- 
tively, and the redshift is z. Emitter's and observer's frame 
of reference are indicated by subscripts 'em' and 'obs'. The 
value of g signals redshift (g < 1), blueshift (g > 1), and 
no net shift (g = 1). We characterize simulated line profiles 
by the line core energy (or line centroid) 

.ycorc r-^ j-, ) \/-J 

i.e. the g-factor that belongs to the flux weight of the whole 
line. This is a suitable quantity to quantify redshifts and 
to investigate the transition from Einsteinian to Newtonian 
gravity . Gravitational redshift in the weak field regime es- 
tablishes pure shifts of spectral features without changing 
their intrinsical shape. However, gravitational redshift in the 
strong field regime - strong gravity - produces remarkable 
distortions of spectral shapes if the rest frame feature is 
compared to its analogue in the observer's frame. Distortion 
is a key feature of relativistic spectra exhibiting very skewed 
and asymmetric line profiles (Fabian et al. 1989; Popovic 
et al. 1995; Tanaka et al. 1995). These effects are impor- 
tant only very close to the black hole. The emitter model 
consists of stationary Keplerian rotating and geometrically 
thin rings that lie in the black hole's equatorial plane. In the 
course of studying gravitational redshift it is better to avoid 
blueshift effects. Hence, a face-on orientation of the emit- 
ting region, i = 0°, is a reasonable choice. In that way, the 
logitudinal Doppler effect is suppressed because there is no 
relative emitter motion along the line of sight. However, we 
still have to take into account the transverse Doppler effect 
in this case. This second order GR effect comes into play 
for relativistic speeds of the emitter. We set an inclination 
angle of i = 1° for numerical reasons. Further, other pa- 
rameters are chosen to be black hole spin at Thome's limit 
a = 0.998 M (Thorne 1974) and a prograde (+) Keplerian 
velocity field of emitters, = y/M /(Vr* + ay/M). The 
emissivity for the rings is modelled by a Gaussian radial 
emissivity profile that is endowed with an emission peaking 
at i?pcak- It is e(r) oc exp (-(r - R pca k) 2 /^) as intro- 
duced by MC04. The parameter er r controls the width of the 
Gaussian or the size of the emitting region and is chosen to 
be 0.2 for all runs. A localized Gaussian emissivity set in 
this way mimics a thin and narrow luminous annulus with 
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Fig. 1 Decay of gravitational redshift z with increasing distance 
to the rotating black hole with a — 0.998 M. The ray tracing sim- 
ulations (filled circles) are compared to the analytical redshifts in 
Schwarzschild and the Kerr geometry, respectively (see text for 
analytical expressions). The inclination angles of the rings amount 
to 1°. The error bars reflect that ray tracing delivers slightly vary- 
ing results depending on numerical resolutions for emitting area 
and/or spectrum. 



~ 1 r g distance between inner and outer edge. In our pro- 
cedure the rings are shifted from large to small distances to 
the black hole. We assume that the radial range of interest 



is 2 



< R 



peak 



< 100000 r„. For each simulated emit- 



ting ring we determine the line core redshift z coro by using 
Eqs. (2) and (1). If we assign i? pea k to each value of z cor e, 
we are able to quantify the gravitational redshift as a func- 
tion of distance to the rotating black hole, i.e. the redshift 
gradient. 



3 Two modes of gravitational redshift: shift 
and distortion 

Fig. 1 displays the core redshifts of emission lines, z coro , as 
a function of i? pC ak for rings inclined to i = 1°. The red- 
shifted lines are dominated by gravitational redshift by con- 
struction. The curve shows that the redshift approaches z — > 
0, i.e. that g —> 1, at distances of a few thousand gravita- 
tional radii from the black hole. This is the regime of nearly 
flat space-time and Newtonian physics. But approaching 
the black hole, space-time curvature becomes more signif- 
icant: z grows rapidly and g approaches zero. We compare 
the simulation results with analytical expressions for the 
redshift z. These redshifts follow from Eq. (1) and the fact 
that gi=o° = a/T for a face-on situation, i = 0° (c.f. Lind 
& Blandford 1985). Here a denotes the lapse function of 
the Kerr geometry and T is the Lorentz factor. The redshift 
dependence on T accounts for the transverse Doppler effect 
that is important even for rings tilted to face-on orientation. 
The lapse function a in Boyer-Lindquist form holds (if re- 
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stricted to the equatorial plane, 9 = tt/2; see MC04, Eq. 
3): 

j r{r 2 -2Mr + a 2 ) 
a -\j r 3 +a 2 (r + 2M) ' (3) 

This function reduces to the famous Schwarzschild factor 
a = d\ — 2M/r for a = 0. The emitter has a simple ve- 
locity field represented by a Keplerian rotation so that the 
Lorentz factor reads 

r= (i-„(*> v (*>)~ 1/2 . (4) 

The orbital velocity is expressed in the locally non-rotating 
frame (Bardeen 1970a, 1970b) by 

« W = f(^-w), (5) 

with cylindrical radius u) and frame-dragging frequency w. 
Finally, we therefore arrive at redshift Zi = o° ( r ) = (T — 
a) /a that can be expanded into a binomial for r S> GM/ c 2 . 
Then, we get z oc r -1 in the far field. This is true for both, 
Kerr and Schwarzschild black holes. The resulting redshift 
curve for a non-rotating black hole with a — (dashed 
curve) versus a rotating black hole with a = 0.998 M (dot- 
ted curve) in Fig. 1 show that deviations due to frame- 
dragging become necessary only very close to the black 
hole 1 . We return to this point in the next section. The ray 
tracing simulations confirm the expectations based on sim- 
ple analytics: The gravitational redshift decays with increas- 
ing distance. But all classical black hole solutions are asymp- 
totically fiat. Hence, there is no finite distance at which the 
gravitational redshift vanishes exactly. However, observabil- 
ity poses a limit to what is practical since the resolving 
power of a spectrograph constrains the amount of gravi- 
tational redshift that can be detected. Assuming a spectral 
resolution of 0. 1 A for H/3 as is obtainable, e.g. by VLT in- 
struments, the corresponding critical value of the g-factor is 
g = 0.999979. In our model this shift occurs at a radius of ~ 
75000 r g . For supermassive black holes (SMBHs) of 10 7 - 
10 8 M Q this radius corresponds to 0.05-0.5 pc, whereas for 
stellar-mass black holes of ~ 10 M Q it is 0.01 AU. Though, 
an observational complication is that at large distances to 
the black hole the line shifts are not dominated by gravi- 
tational redshift. Gravitational redshift not only occurs as 
a pure shift: relativistic emission lines are also deformed in 
intrinsic shape. But this mode only occurs in sufficient prox- 
imity of the emitter to the black hole. Hence, we refer to this 
as strong gravity. It is suitable to change the inclinations to 
higher angles to demonstrate this second effect. Therefore, 
we set i = 75° (any other simulation parameters remain) 
and investigate the characteristic broad line profile with the 
two relic Doppler peaks for various distances of the rings 
to the black hole. We make use of two quantities to illus- 
trate the distortion effect. The first quantity termed g lp is 
the (/-factor associated with the red relic Doppler peak that 
is shifted to lower energies as the ring approaches the hole, 

1 In this Keplerian model plausible redshifts are only computed for radii 
greater than the marginally stable orbit, r ms (a) 
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Fig. 2 Distortion of the red relic Doppler peak for rings satisfy- 
ing i = 75°: With decreasing radius the red peak can be found at 
lower peak energies due to strong gravitational redshift. Very close 
to the black hole - r < 3 r g at this specific inclination - these dis- 
tortions are such strong that the red Doppler peak is highly blurred 
and vanishes effectively in the line profile. 
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Fig. 3 Doppler peak spacing (DPS) in units of g derived from 
line emitting rings at i — 75° . Far away from the black hole both 
peaks approach significantly and the double-peaked line profile 
becomes very narrow: At ~ 3000 r g the peak difference in g only 
amounts to ~ 0.03. 



displayed in Fig. 2. This plot can be used to read the half- 
energy radius (i.e. the radius where the g-factor associated 
with the red Doppler peak is exactly 1/2) associated with 
the red relic Doppler peak. At i = 75° this quantity can be 
found at ~ 5 r g and documents that strong gravity is im- 
portant only very close to the black hole. The distortion of 
the intrinsic line profile by gravitational redshift can be seen 
from g Tp alone or from the second quantity called Doppler 
peak spacing (DPS) that is computed from the g-factors 
associated with the red peak (rp) and the blue peak (bp): 
DPS = gb P — 5r P - Fig. 3 shows that the spacing between 
the Doppler peaks does not remain constant as the emitter 
approaches the black hole: The DPS value rises quickly, i.e. 
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the line is stretched. Gravitational redshift causes an addi- 
tional suppression in flux so that close to the black hole the 
emission line from an intermediately to highly inclined ring 
is asymmetric and skewed. 

4 Effects of black hole rotation 

A rotating black hole gains support from black hole growth 
theories: Shapiro (2005) reviewed spin values of SMBHs 
in the local universe to be at either a ~ 0.95 M (MHD 
disk) or even a ~ M (standard thin gas disk). But rotation 
of space-time decreases very steeply with distance to the 
black hole, to oc r~ 3 . Fig. 1 already documents a difference 
between the Kerr and the Schwarzschild redshifts only for, 
r < 10 r g . Therefore, black hole rotation can only be probed 
with spectral features originating in regions very close to the 
black hole, like X-ray iron K lines. In contrast, optical BLR 
emission lines in AGN are not suited for probing black hole 
rotation in that manner. These optical and X-ray lines to- 
gether cover a radial scale of five orders of magnitude. Only 
larger multi-wavelength data sets for each particular source 
allow for testing gravitationally redshifted lines in differ- 
ent spectral ranges. All these data should point towards the 
same central mass that could be determined this way. High- 
quality data from the black hole's vicinity could even give 
rise for its spin. Kollatschny (2003) and also Miiller & Wold 
(2006) demonstrated for the Narrow-line Seyfert-1 galaxy 
Mrk 110 that gravitational redshift studies also reveal the 
inclination angle of the inner accretion disk if the flat, Kep- 
lerian BLR scenario holds. This offers a new powerful black 
hole mass estimator besides reverberation mapping, M-a 
relation, stellar velocity dispersion or Maser studies. 

5 Conclusions 

We analyse emission line cores that originate from Keple- 
rian rotating rings at 2 to 100000 r g distance to a rapidly 
spinning black hole by using relativistic Kerr ray tracing 
simulations. These simulations confirm the expectation that 
the redshift gradient follows a simple z w ,. , 1 — n law as 

° r distancc[r g J 

a straightforward consequence of the Schwarzschild factor. 
However, closer to the black hole the behaviour is not that 
simple and redshift is influenced by non-linear effects and 
black hole spin. Generally, gravitational redshift occurs in 
two modes: One regime starts at larger distances from the 
black hole and shifts only the line as a total feature while 
conserving its intrinsic shape. The other regime, which is 
the strong gravity regime, dominates at r < 10 r g . Emis- 
sion lines that form in the vicinity of the black hole are 
strongly deformed and suppressed. GR predicts that grav- 
itational redshift is a ubiquitous phenomenon. However, the 
underlying redshift may be blurred by blueshift effects, e.g. 
from outflows. If the amount of blueshift is known or negli- 
gible, the verification of gravitational redshift is only detector- 
dependent. Assuming currently feasible optical resolving 



powers, we find a limiting distance of 75000 r g where grav- 
itationally redshifted optical lines can be probed. But com- 
peting effects and more complex physics are likely to be in- 
volved so that the verification remains challenging. We sug- 
gest to intensify high-resolution multi-wavelength obser- 
vations in varying distance around black hole candidates so 
that the central mass can be derived from systematic spectral 
shifts. If our flat Keplerian model applies for line emitting 
regions, it is even possible to determine the inner inclination 
angle. Black hole spin can only be probed by spectral fea- 
tures that form sufficiently close to the center, e.g. by iron 
Ka lines. Our study holds for any black hole of arbitrary 
mass, i.e. stellar-mass black holes in XRBs, IMBHs that 
may be found in ultra-luminous X-ray sources or globular 
clusters, and SMBHs in the cores of galaxies. 
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